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Abstract 

In addition to three active neutrinos v e , v„ and v T , one or more light sterile 
neutrinos have been conjectured to account for the LSND, MiniBooNE and 
reactor antineutrino anomalies (at the sub-eV mass scale) or for warm dark 
matter in the Universe (at the keV mass scale). Heavy Majorana neutrinos 
at or above the TeV scale have also been assumed in some seesaw models. 
Such hypothetical particles can weakly mix with active neutrinos, and thus 
their existence can be detected at low energies. In the (3+3) scenario with 
three sterile neutrinos we present a full parametrization of the 6x6 flavor 
mixing matrix in terms of fifteen rotation angles and fifteen phase angles. 
We show that this standard parametrization allows us to clearly describe 
the salient features of some problems in neutrino phenomenology, such as 
(a) possible contributions of light sterile neutrinos to the tritium beta decay 
and neutrinoless double-beta decay; (b) leptonic CP violation and deformed 
unitarity triangles of the 3x3 flavor mixing matrix of three active neutrinos; 
(c) a reconstruction of the 6x6 neutrino mass matrix in the type-(I+II) 
seesaw mechanism; and (d) flavored and unflavored leptogenesis scenarios in 
the type-I seesaw mechanism with three heavy Majorana neutrinos. 
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I. INTRODUCTION 



One of the fundamental questions in neutrino physics and cosmology is whether there 
exist extra species of neutrinos which do not directly participate in the standard weak 
interactions. Such sterile neutrinos are certainly hypothetical, but their possible existence 
is either theoretically motivated or experimentally implied. For example, 

• heavy Majorana neutrinos at or above the TeV scale are expected in many seesaw 
models [1], which can not only interpret the small masses of three active neutrinos but 
also account for the cosmological matter-antimatter asymmetry via the leptogenesis 
mechanism [2]; 

• the LSND antineutrino anomaly [3] , the MiniBooNE antineutrino anomaly [4] and the 
reactor antineutrino anomaly [5] can all be explained as the active-sterile antineutrino 
oscillations in the assumption of two species of sterile antineutrinos whose masses are 
close to 1 eV [6]; 

• an analysis of the existing data on the cosmic microwave background (CMB), galaxy 
clustering and supernovae la favors some extra radiation content in the Universe and 
one or two species of sterile neutrinos at the sub-eV mass scale [7] 1 ; 

• sufficiently long-lived sterile neutrinos in the keV mass range can serve for a good 
candidate for warm dark matter, whose presence may allow us to solve or soften 
several problems that we have recently encountered in the dark matter simulations [9] 
(e.g., to damp the inhomogeneities on small scales by reducing the number of dwarf 
galaxies or to smooth the cusps in the dark matter halos) 2 . 

No matter how small or how large the mass scale of sterile neutrinos is, they are undetectable 
unless they mix with three active neutrinos to some extent. The strength of active-sterile 
neutrino mixing can be described in terms of some rotation angles and phase angles, just 
like the parametrization of the 3x3 quark flavor mixing in the standard model [14]. 

The main purpose of this paper is to present a full parametrization of the 6x6 flavor 
mixing matrix IA in the (3+3) scenario with three sterile neutrinos denoted as u x , v and v z : 
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x If the bound obtained from the Big Bang nucleosynthesis is taken into account, however, only 
one species of light sterile neutrinos and antineutrinos is allowed [8] . 

2 There are some interesting models which can accommodate sterile neutrinos at either keV [10] 
or sub-eV [11] mass scales. A model-independent argument is also supporting the conjecture of 
warm dark matter particles hiding out in the "flavor desert" of the fermion mass spectrum [12]. 
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where v i (for % = l,---,6) stand for the mass eigenstates of active and sterile neutrinos. 
Such a complete parametrization, which has been lacking in the literature [13], is expected 
to be very useful for the study of neutrino phenomenology at both low and high energy 
scales. We propose a simple but novel way to establish the connection between active and 
sterile neutrinos in terms of fifteen mixing angles and fifteen CP-violating phases. It allows 
us to clearly describe the salient features of some interesting problems, such as (a) possible 
contributions of light sterile neutrinos to the tritium beta (/?) decay and neutrinoless double- 
beta (Oz/2/3) decay; (b) leptonic CP violation and deformed unitarity triangles of the 3x3 
flavor mixing matrix of three active neutrinos; (c) a reconstruction of the 6x6 neutrino mass 
matrix in the type-(I+II) seesaw mechanism; and (d) flavored and unflavored leptogenesis 
scenarios in the type-I seesaw mechanism with three heavy Majorana neutrinos. 



II. THE STANDARD PARAMETRIZATION 



The 6x6 unitary matrix U defined in Eq. (1) can be decomposed as 



U 
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(2) 



in which and 1 stand respectively for the 3x3 zero and identity matrices, U and V are 
the 3x3 unitary matrices, and A, B, R and S are the 3x3 matrices which satisfy the 
conditions 



AA^ + RR^ 
AS^ + RB^ 
A^A + S^S ■ 



■ BB^ + SS^ 
A^R + S ] B 
B^B + FtfR 



1, 
0, 
1 , 



(3) 



as a result of the unitarity of U. In the limit of R = S = 0, A — B — 1 holds and thus 
there is no correlation between the active sector (described by V ) and the sterile sector 
(characterized by U ). In view of Eq. (A7) in Appendix A, we parametrize U as follows: 
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where fifteen two-dimensional rotation matrices (for 1 < i < j < 6) in a six-dimensional 
complex space have been given in Eqs. (A2) — (A6). To be explicit, 
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in which ■ = cos 9^ and = e l5 ^ sin 6^ with ^ ■ and 8^ being the rotation angle and phase 
angle, respectively. Both V and U have the standard form as advocated in Ref. [14], and 
either of them consists of three mixing angles and three CP-violating phases. If the sterile 
sector is switched off, we are then left with the 3x3 unitary matrix V which describes the 
flavor mixing of three active neutrinos. If the active sector is switched off, one will arrive 
at the 3x3 unitary matrix U which purely describes the flavor mixing of three sterile 
neutrinos. In the type-I seesaw mechanism [1], for example, U is essentially equivalent to 
the unitary transformation used to diagonalize the 3x3 heavy Majorana neutrino mass 
matrix M R and therefore relevant to the leptogenesis mechanism [2]. 

With the help of Eq. (4) and Eqs. (A2) — (A6), a lengthy but straightforward calculation 
leads us to the explicit expressions of A, B, R and S as follows: 
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V +Si4S 1 5S2 5 S2 6 C 36 — Si4C 15 S 16 C 26 C 36 
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We see that the textures of A and B are rather similar, so are the textures of R and S. In 
fact, the expression of B can be obtained from that of A* with the subscript replacements 
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15 -H- 24, 16 -H- 34, and 26 <r> 35; and the expression of S can be obtained from that of — R* 
with the same subscript replacements. Note that the results of A and R have been obtained 
in Ref. [15], and here we present the results of B and S to complete a full parametrization 
of the 6x6 flavor mixing matrix U. 

It proves convenient to define V = AV and U = C/ S which describe the flavor mixing 
phenomena of three active neutrinos and three sterile neutrinos, respectively. Furthermore, 
S = C/o^^o links the mass eigenstates (z/ 1? u 2 , u 3 ) to the sterile flavor eigenstates (v x , v , u z ) 
in the chosen basis. We therefore have 



V 



Vn 



+ R 



(8) 



and 



(9) 

Eq. (8) directly leads us to the standard weak charged- current interactions of six neutrinos: 

(10) 
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where V is just the Maki-Nakagawa-Sakata-Pontecorvo (MNSP) matrix [16] responsible for 
the active neutrino mixing, and R measures the strength of charged-current interactions 
between (e, //, r) and (u 4 ,v 5 ,u 6 ). Because of 

VV ] = AA^ = 1 - RB* , 

rfU = B ] B = 1- R ] R, (11) 

we find that both V and U are not exactly unitary and their non-unitary effects are simply 
characterized by non-vanishing R and S. 

In view of current observational constraints on sterile neutrinos, we expect that the 
mixing angles between active and sterile neutrinos are strongly suppressed (at most at the 
0(0.1) level [6] 3 ). The smallness of 9^ (for i — 1, 2, 3 and j = 4, 5, 6) allows us to make the 
following excellent approximations to Eqs. (7) and (8): 
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3 For example, the non-unitarity of V = AV or the deviation of V from V can at most be at the 
1% level as constrained by current neutrino oscillation data and precision electroweak data [17]. 
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where the terms of 0(sf ) have been omitted; and 
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where the terms of 0(s?) have been omitted. It turns out that R ~ — holds in the same 
approximation. 

Note that the 6x6 unitary matrix U can be used to describe not only the flavor mixing 
between active and sterile neutrinos but also the flavor mixing between ordinary and extra 
quarks. Note also that it is straightforward to obtain the (3+1) flavor mixing scenario from 
Eqs. (6) and (7) by switching off the mixing angles 9 i5 and # - 6 (for 1 < i < 4 and 1 < j ' < 5), 
or the (3+2) flavor mixing scenario by turning off the mixing angles - 6 (for 1 < j < 5). 



III. SOME APPLICATIONS 



To illustrate the usefulness of our parametrization of the 6x6 flavor mixing matrix U, 
let us briefly discuss its four simple but instructive applications in neutrino phenomenology. 



A. The effective masses of j3 and 0v2/3 decays 



One or two light sterile neutrinos at the sub-eV mass scale have been hypothesized 
for a quite long time to interpret the LSND antineutrino anomaly [3] and the subsequent 
MiniBooNE antineutrino puzzle [4]. In general, however, it seems more natural to assume 
the number of sterile neutrino species to be equal to that of active neutrino species [18] 4 , 
such that even possible warm dark matter in the form of one or two species of keV sterile 
neutrinos could be taken into account. 

For simplicity and illustration, we are only concerned about the effective masses of the 
tritium beta ((3) decay — > pie + e~ + V e and the neutrinoless double-beta (Oz/2/3) decay 
A(Z, A) — > A(Z + 2, iV — 2) + 2e~ in the (3+3) neutrino mixing scenario. The former is 



(m)' e 



1/2 



— \[{ m )i c li c lb c i& + m 4 s 14 c 15 c 16 + m 5 s 15 c 16 + m 6 s 16 



\2^2 
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(14) 



where (m) e = ^micl 2 cl 3 + ^2 s i2 c i3 ~l~ m 3 s i3 is t-h- e standard contribution from three active 
neutrinos. We see that (m)' e > (m) e always holds. The effective mass of the 0u2/3 decay is 



4 In order to avoid any severe conflict between such a (3+3) scenario and the standard ACDM 
cosmology, it is perhaps necessary to either loosen the mass hierarchy of three sterile neutrinos 
(i.e., not all of them are of O(0. 1) eV) or refer to some nonstandard models of cosmology [7]. 
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( m )'ee = H m i V ei = ( m )ee (c 14 C 15 C 16 f + m 4 (s* l4 C lb C 1Q f + m 5 (st 5 C 16 ) 2 + m & (s* 6 ) 2 (15) 
i=l 

with (m) ee = ^i(c 12 c 13 ) 2 + ^2(^i2 c i3) 2 + m 3(^h) 2 being the standard contribution from 
three active neutrinos. It is difficult to say about the relative magnitudes of (m) ee and (m)' ee , 
because the CP-violating phases may give rise to more or less cancelations of different terms 
in them. In particular, even (m) ee = [19] or (m)' ee = [20] is not impossible. If both 
(m)' e and (m)' ee can be determined or constrained in the future experiments, a comparison 
between them might be able to probe the existence of light sterile neutrinos [21]. 



B. Deformed unitarity triangles and CP violation 

Switching off three sterile neutrinos, one may describe flavor mixing and CP violation of 
three active neutrinos in terms of six unitarity triangles in the complex plane [22]. Three of 
them, defined by the orthogonality conditions 



A, 



Wi + v% + v% = > 

KiKi + v r2 v; 2 + v r3 v; 3 = o , 
v A v^ + v e2 v; 2 + v e3 v; 3 = o , 



(16) 



are illustrated in FIG. 1 (left panel). The area of each triangle is equal to J /2, where J is 
the Jarlskog parameter given in Eq. (B2) and measures the strength of leptonic CP-violating 
effects in v — > is T , v T — > v e and v e — > v oscillations. Now let us turn on the contributions 
of three sterile neutrinos to flavor mixing and CP violation. Then Eq. (16) approximates to 



a; 
a; 

A' 



^iK*i + v, 2 v; 2 + v^ 3 v; 3 ~ -z* 
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V el V* + W + V e3 V L 



jj,3 



-X* 



(17) 



where X = s 14 s* 24 +s 15 s* 25 +s 16 s* 26 , y = hi s li + hi s l5 + " s i^m and Z = s 24 s* 34 +s 25 s* 35 +s 26 s* m . 
These deformed unitarity triangles are also illustrated in FIG. 1 (right panel). The small 
differences of their areas from J /2 just signify the new CP- violating effects. 

Let us take a look at the CP-violating asymmetries between u a — > and V a — > V ^ 
oscillations, defined as A a/3 = P{y a — > v^) — P(v a — > V^). With the help of Eq. (B5), we 
explicitly obtain 



\e ~ - 4 ( J + C 12 S 12 C 23 lmX) sin , 



At - - 4 ( J o + c 12 s 12 s 23 lmY) sin 



2E 
Am^L 

2E 



\r - +4 [Jq + C 12 S 12 C 23 S 23 (s 23 ImX + C 23 ImF )] Sm 



Am 21 L 



2E 

-i<5 



+ 4c 23 s 23 lmZ sin 



Am\ 2 L 
2E 



, (18) 



where X = Xe l5 ^ i y = y e ^(s 12 +s 23 ) and Z = Ze td ^. We see that A^ e and A er are 
related to the deformed unitarity triangles A' T and A' respectively. In comparison, A^ 
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has something to do with A„ T . It is therefore possible to determine three new CP- violating 
terms ImX, ImY and \m.Z by measuring the CP- violating effects in neutrino oscillations. 
Note that three CP-violating asymmetries in Eq. (18) satisfy the correlation 

\r + {4s\e + 4sAr) ~ ^S^lmZ sin . (19) 

When Aml 2 L/E ~ 7r holds, both A^ e and A er are suppressed such that A^ T becomes a 
pure measure of the non-unitary CP- violating parameter ImZ. This interesting possibility, 
together with terrestrial matter effects, has been discussed before (e.g., Refs. [15] and [26]). 



C. Reconstruction of the 6x6 neutrino mass matrix 

The type-(I+II) seesaw mechanism [27] is a good example to illustrate the flavor mixing 
between three active neutrinos and three heavy Majorana neutrinos. In this mechanism the 
mass term of six neutrinos is usually written as 

-£-=5 J£ )(£) + !«., (20) 

where u L and N R represent the column vectors of three left-handed neutrinos and three 
right-handed neutrinos, respectively. The overall 6x6 neutrino mass matrix in Eq. (20) 
can be diagonalized by a unitary transformation: 

where U is already given in Eq. (2), M v = Diag{m 1 , m 2 , fn 3 } and M N = Diag{M 1 , M 2 , M 3 } 
with m i or M i (for % — 1, 2, 3) being the physical masses of light or heavy Majorana neutrinos. 
The standard weak charged-current interactions of six neutrinos are given by Eq. (10) with 
z/ 4 = M 1 , z/ 5 = N 2 and v 6 = N 3 in the basis of mass eigenstates. With the help of Eq. (2), 

M L = VM V V T + RM N R T ~ V Q M v Vg + RM N R T , 
M D = VMjF + RM N U T ~ RM n Uq , 

M R = SM V S T + UM N U T ~ U M N U£ , (22) 

where V = ^V'q, [/ = C/ S and 5 = C/qS'V^ have been defined before. The approximations 
made on the right-hand side of Eq. (22) follow the spirit that only the leading terms of M L , 
M D and M R are kept. It is then possible to reconstruct these 3x3 neutrino mass matrices, 
at least in principle, in terms of neutrino masses and flavor mixing parameters [28]. 

Given the basis where M R is diagonal, real and positive, Eq. (22) implies that M R ~ M N 
together with U ~ 1 is a good approximation. Note that such a flavor basis is often chosen 
in the study of leptogenesis, because the decays of N i (for % = 1,2,3) need to be calculated. 
It is also much easier to reconstruct M D and M L in this special basis. For example, 



M D ~ RM 



N 



(M^U M 2 s* 15 M 3 s* 16 

M ± s* 24 M 2 s* 25 M 3 s* 26 | ; (23) 
{MA M 2 s* 35 M 3 s* 6 , 
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and six independent matrix elements of M L ~ V M U V^ + M D R T can similarly be obtained: 

(M L ) ee ~ m x (c 12 c 13 ) 2 + m 2 (s* l2 c l3 ) 2 + m 3 (s* 3 ) 2 + M x {s\ 4 f + M 2 (s* l5 ) 2 + M 3 (s* 6 ) 2 , 

\M\)eti — ~ m l C 12 C 13 (^12 C 23 + C 12^13^23) + m 2^12 C 13 ( C 12 C 23 — ^12^13^23) + m 3 C 13^13^23 

+M 1 s* 14 s 24 + M 2 s* 15 s* 25 + M 3 s* 16 s 26 , 

(MiJ eT — m \ c l2 C 13 (^12^23 — C 12^13 C 23) — m 2^12 C 13 ( C 12^23 + ^12^13 C 23) + m 3 C 13^13 C 23 

+M 1 s* 14 s* 34 + M 2 s* 15 s* 35 + M 3 s* 16 s* m , 
— tti 1 (s 12 c 23 + c 12 s 13 s 23 ) + m 2 (c 12 c 23 — s 12 s 13 s 23 ) + m 3 (c 13 s 23 ) 
+M X (s* 24 ) 2 + M 2 (s* 25 ) 2 + M 3 (s* 26 ) 2 , 
{Ml)^ t — —m 1 (s 12 c 23 + c 12 s 13 s 23 ) (s 12 s 23 — c 12 s 13 c 23 ) 

— m 2 (,C 12 C 2 3 — S 12 S 13 S 2 3J l C 12 5 23 5 12 5 13 C 23/ m 3 C 13 C 23 5 23 

+M 1 s* 4 s* 34 + M 2 s* 2b s* 35 + M 3 s* 26 s* 36 , 

(-^l)tt — m l (^12^23 — C 12^13 C 23) + m 2 ( C 12^23 + ^12^13 C 23) + m 3 ( C 13 C 23) 

+M 1 (s* 34 ) 2 + M 2 (s* 5 ) 2 + M 3 (s* 36 ) 2 . (24) 

So any specific textures of M D and M L predicted in a specific type-(I+II) seesaw model must 
have direct and important impacts on the magnitudes of neutrino masses, flavor mixing 
angles and CP-violating phases. 



D. Flavored and unflavored leptogenesis scenarios 

It is straightforward to obtain the type-I seesaw mechanism from the type-(I+II) seesaw 
mechanism by switching off the M L term. In this special case one arrives at the exact type-I 
seesaw relation VM U V T + RM N R T = 0, where V and R satisfy the unitarity condition 
VV^ + RR^ = 1. Given M i 3> m-, it is more popular to write the 3x3 light Majorana 
neutrino mass matrix as 



VM V V 2 



-RM N R 



-m d m^mE , 



(25) 



where V = AV ~ V holds in the same approximation [29]. Associated with this seesaw 
picture, the leptogenesis mechanism [2] provides a natural possibility of accounting for the 
observed matter-antimatter asymmetry of the Universe via the lepton-number-violating, 
CP-violating and out-of-equilibrium decays of N i and the (B — L)-conserving sphaleron 
processes. The CP-violating asymmetry between — >■ £ a + H and N { — >■ I a + H decays, 
denoted as e ia (for i — 1,2,3 and a — e, /i,t), is given by [30] 



TjNj ^i a + H)- r(iV t ^£ a + H) 

r(7V i ^4 + //) + r(7V.^4 + ^) 



8tto 2 (MjM D ; 



+ Im 



(M*) Qi (M D ) QJ .(MtM D )* J g 




(26) 



9 



where the loop functions J-(x) = \/x{(2 — x)/(l — x) + (1 + x) \a[x/ (1 + x)]} and Q{x) = 
1/(1 — x) have been introduced. If all the interactions in the era of leptogenesis were blind 
to lepton flavors, then only the total CP- violating asymmetry e i should be relevant: 



i 



8tt(M£M d ) 



E Im [(j^d) 




(27) 



The leptogenesis mechanisms associated with Eqs. (26) and (27) are usually referred to as 
flavored and unflavored leptogenesis scenarios, respectively. 

In the flavor basis where Mr ~ M N and U ~ 1 hold, we have obtained the explicit 
expression of M D ~ RM N in Eq. (23). It is then straightforward to arrive at 



(Mi M D ) n ~ M 1 2 ( J Rt jR)n ~ Ml (4 + 4 + 4 
(M*M D ) 22 ~ Mi(rfR) 22 ~ M 2 2 (s 2 5 + 4 + 4 
(Mt M D ) 33 ~ M 3 2 ( J Rt jR ) 33 ~ M 3 2 (4 + 4 + S 2 



2 ") ■ 

36 y ) 



(28) 



and 



'15 



'24°25 



'34°357 ) 



(MjM D ) 12 ~ M 1 M 2 (i?t jR ) 12 ~ Ml M 2 (s 14 s; 

(mJm d ) 13 ~ M 1 M 3 (i?t jR)i3 ~ Mi m 3 (s 14 s; 6 + 44 + 44) , 

(MiM D ) 23 ~ M 2 M 3 (i? t J R) 23 ~ M 2 M 3 (s 15 ^ 6 + 44 + 44) . 



(29) 



So the unflavored CP-violating asymmetry e i depends on nine phase differences 8 iA — 5, 



5 i4 — 5 i6 and 5 i5 — 5 i6 (for i = 1,2,3), but only six of them are independent. Because V ~ V 
holds in the same approximations as made above, we conclude that there is not any direct 
relationship between the CP violation at low energies (governed by 5 12 , 5 13 and 5 23 ) and 
the unflavored leptogenesis at high energies. As for the flavored leptogenesis, the relevant 
CP- violating asymmetries e ia also rely on the aforementioned nine phase differences. This 
point can be clearly seen from 



Im 
Im 



(M*) Qi (M D ) Qj (MiM D )J ~ M 2 M 2 Im [W ai R aj {I$ R) i3 
(M*) m (M D 4(Mt M D yJ ~ M 2 M 2 Im [R* ai R a .(tfR)T. 



(30) 



Given the masses of three heavy Majorana neutrinos and the flavor mixing parameters 
between light and heavy Majorana neutrinos, it is then possible to determine e ia and e i (for 
2 = 1,2,3 and a — e, /i, r) so as to realize the flavored or unflavored leptogenesis mechanism. 



IV. SUMMARY AND REMARKS 

An appealing and puzzling feature of the standard model is that it happens to have three 
species of leptons and quarks. If extra species of matter particles exist, no matter whether 
they are sequential in or exotic to the standard model itself, they definitely signify new 
physics. In this paper we have conjectured the presence of three species of sterile neutrinos 
and presented a full parametrization of the 6x6 flavor mixing matrix for active and sterile 
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neutrinos in terms of fifteen rotation angles and fifteen phase angles. Such an exercise makes 
sense because we do have some preliminary observational hints on light sterile neutrinos, and 
the theoretical motivation for the existence of heavy Majorana neutrinos in the seesaw and 
leptogenesis mechanisms is also very strong. 

We have shown that this standard parametrization of the 6x6 flavor mixing matrix in 
the (3+3) scenario of active and sterile neutrino mixing allows us to clearly describe the 
salient features of some problems in neutrino phenomenology. Four examples have been 
briefly discussed in this connection: (a) possible contributions of light sterile neutrinos to 
the tritium beta decay and neutrinoless double-beta decay; (b) leptonic CP violation and 
deformed unitarity triangles of the 3x3 flavor mixing matrix V of three active neutrinos; 
(c) a reconstruction of the 6x6 neutrino mass matrix in the type-(I+II) seesaw mechanism; 
and (d) flavored and unflavored leptogenesis scenarios in the type-I seesaw mechanism with 
three heavy Majorana neutrinos. Our results are expected to be useful to understand the 
impacts of extra neutrino species on the standard weak interactions and neutrino oscillations 
in a better and more straightforward way. Furthermore, the parametrization itself can also 
be applied to the quark sector if extra species of quarks are conjectured. 

Let us stress that the presence of new degrees of freedom in the neutrino sector may 
apparently violate the unitarity of the 3x3 flavor mixing matrix V of three active neutrinos 
in the weak charged-current interactions. Hence testing the unitarity of V in neutrino 
oscillations and searching for the signatures of new neutrinos at TeV-scale colliders can 
be complementary to each other, both qualitatively and quantitatively, in order to deeply 
understand the intrinsic properties of Majorana particles. On the other hand, light or heavy 
sterile neutrinos can have significant consequences in cosmology. The latter provides us 
with a good playground to study hot dark matter in the presence of sub-eV sterile neutrinos 
and warm dark matter in the form of keV sterile neutrinos, together with the cosmological 
matter-antimatter asymmetry via leptogenesis in which heavy Majorana neutrinos and their 
decays play the key role. We hope that any experimental breakthrough in these aspects will 
pave the way towards the true theory of massive neutrinos. 

The author would like to thank B. Dziewit, Y.F. Li, W. Rodejohann, S. Su, S. Zhou and 
M. Zralek for useful discussions. This work was supported in part by the National Natural 
Science Foundation of China under grant No. 11135009. 
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APPENDIX A 



The 6x6 unitary matrix U in Eq. (1) can be expressed as a product of fifteen two- 
dimensional rotation matrices in a six-dimensional complex space [15]: 

U = (0 56 46 36 26 16 ) (0 A5 35 25 O l5 ) (0 34 24 14 ) (0 23 13 ) 12 , (Af ) 

where O tj (for 1 < i < j < 6) are unitary and read as follows: 
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In the above equations we have defined q ■ = cos 9^ and = e % a sin 0^ with Q { - and 8^ 
being the rotation angle and phase angle, respectively. Because some of the two-dimensional 
rotation matrices can commute with each other, it is possible to rewrite Eq. (Al) as 

U = (0 56 46 45 ) (0 36 26 16 35 25 15 O u 2A O u ) (0 23 13 12 ) . (A7) 

On the right-hand side of Eq. (A7), the combination {0 23 13 O l2 ) purely describes the flavor 
mixing among three active neutrinos while (0 56 46 45 ) purely describes the flavor mixing 
among three sterile neutrinos. It is the combination (0 36 ■ • ■ 14 ) in Eq. (A7) that allows 
the active and sterile sectors to "talk" to each other, as discussed in section II. 



APPENDIX B 

Given the n x n flavor mixing matrix, one may always calculate its (n — l) 2 (n — 2) 2 /4 
rephasing invariants of CP violation, the so-called Jarlskog parameters [23] . As for the 6x6 
unitary matrix U under discussion, we totally have 100 invariants of this nature. But we 
are mainly concerned about the Jarlskog parameters of V = AV defined in section II: 

J% = Im(V^.V^) , (Bl) 

in which the Greek indices run over (e,/i, r) and the Latin indices run over (1,2,3). Since 
V describes the flavor mixing of three active neutrinos, J^g measure the corresponding CP- 
violating effects in their oscillations. In the absence of three sterile neutrinos (i.e., A = 1 
and V = V ), one arrives at a universal Jarlskog parameter 

^0 = Jen = Je» = = ^ \xt = ^Jr = ^ftr = ^re = ^re = ^re = C 12 S 12 C 13 S 13 C 23 S 23 ^ $ (^2) 
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with 5 = 5 13 — 5 12 — 5 23 , as guaranteed by the unitarity of V . When the contributions of 
three sterile neutrinos are switched on, J£g can be calculated with the help of Eqs. (5) and 
(12) in a good approximation. One finds 

J ll - J o + c 12 s 12 c 23 lmX , 

J re — Jq + c 12 S 12 S 23-'- m ^ ? 

J ir - J + C 12 S 12 C 23 S 23 ( S 23 ImX + C 23 ImF ) > 

- ^0 + C 12 C 23 S 23 ( S 12 S 23 ImX + «12 C 23 Imy + C 12 ImZ) , 

J^l - J o + S12C23S23 (ci 2 s 2 3 lmX + Ci 2 c 23 Imr - s 12 ImZ) , (B3) 

and ~ J3i ~ J23 ~ j3i ~ j q [24] ? w here X = Xe~ iS ™, Y = y e - i( - s ™ +s ^ and Z = Ze~ iS ™ 
with X, y and Z being defined below Eq. (17). Note that we have assumed 8 13 , 9 i4 , 6 i5 
and 8 i6 (for % = 1,2,3) to be small in our calculations, and thus the terms of (9(s 13 |X|), 
0(s 13 |Y|) and (9(s 13 |Z|) together with those higher-order terms have been omitted from the 
above results. The fact that ~ J^J ~ ~ ^ J holds in the above approximation 
is simply because they all involve the smallest matrix element of V (i.e., V e3 ~ s* 3 ) [25]. 

It is well known that the maximal value of J is J™ ax = l/(6\/3) ~ 9.6% [23]. In 
comparison, the magnitudes of ImX, ImY and ImZ are likely to reach the percent level if 
those active-sterile mixing angles are of 0(0.1) and the relevant CP- violating phases are of 
0(1). So the five Jarlskog parameters in Eq. (B3) might deviate from the standard one J 
in a significant way, depending on the constructive or destructive contributions from three 
sterile neutrinos. It is therefore important to observe all the nine Jarlskog parameters in 
neutrino oscillations. Taking account of the non-unitarity of V, one may easily derive the 
probabilities of u a — > u B and V a — > V B oscillations [15,17]: 



E KiW + 2 £ Re (V^V:^) cos Ay - 2 £ J% sin Ay 



i<j i<j 



£ IKJ'I V + 2 E Re (K.i^V^-VyJ cos Ay + 2 £ J^j sin Ay 

*<*.^> = J ,<J (vy .) (vv) " ' (B4) 

where Ay = Am|L/ (2_E) with Amy- = mf — m|, _E being the neutrino beam energy and L 
being the baseline length. As a consequence, 

4* = P (»a - Pfa = ^ g 4 Sin 

~ 4 [jg sin A 21 + ( + jg) sin A 32 ] , (B5) 

where ~ 1 and A 31 ~ A 32 (i.e., Am 31 ~ Am| 2 [14]) have been taken into account. This 
result implies that both and + J^g can in principle be determined if the baseline of 
neutrino oscillations is sufficiently long. Of course, terrestrial matter effects may more or 
less contaminate the genuine CP-violating effects in such long-baseline neutrino oscillation 
experiments and should be properly treated [26]. 
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FIGURES 
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FIG. 1. Left panel: three unitarity triangles of V = AVq in the absence of three sterile 
neutrinos (i.e., A = 1 and V = V ); Right panel: three deformed unitarity triangles of V = AV 
in the presence of three sterile neutrinos, where X, y and Z are defined below Eq. (17). 
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